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Abstract 

Compton scattering on the deuteron is studied in the framework of baryon chiral 
perturbation theory to third order in small momenta, for photon energies of order the 
pion mass. The scattering amplitude is a sum of one- and two-nucleon mechanisms 
with no undetermined parameters. Our results are in good agreement with existing 
experimental data, and a prediction is made for higher-energy data being analyzed 
at SAL. 
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1 Introduction 

Experimental facilities which accurately measure the energy of a photon beam using 
photon tagging have made possible a new generation of experiments which probe the 
low-energy structure of nucleons and nuclei. In particular, photon tagging can be used to 
measure Compton scattering on weakly-bound systems, since it facilitates the separation 
of elastic and inelastic cross sections. At sufficiently low energy u> the spin-averaged 
forward Compton scattering amplitude for any nucleus is, in the nuclear rest frame: 



where e and e' are the polarization vectors of the initial and final-state photons, Ze is 
the total charge of the nucleus, and AM is its total mass. The first term in this series 
is a consequence of gauge invariance, and is the Thomson limit for low-energy scattering 
on a target of mass number A and atomic number Z. The coefficient of the second term 
is the sum of the target electric and magnetic polarizabilities, a and f3, respectively. The 
polarizabilities can be separated by going to backward angles, which probes the difference, 
a — p. The polarizabilities contain information about the electromagnetic structure of 
the target. 

A number of recent experiments have focused on Compton scattering on the proton 
and on the deuteron. These experiments have already resulted in a wealth of new data on 
a p and /3 P , the proton electric and magnetic polarizabilities, and might yield information 
on a n and /3 n , the corresponding neutron quantities. These are basic low-energy hadron 
parameters which reflect the underlying QCD dynamics. In a simple quark model picture 
they contain averaged information about the charge and current distribution produced by 
the quarks inside the nucleon. 

At energies well below the chiral symmetry breaking scale, A x ~ 4nf w ~ M ~ m p , the 
electromagnetic interactions of pions and nucleons can be described systematically using 
an effective field theory. This effective field theory, known as chiral perturbation theory 
(xPT), reflects the observed QCD pattern of symmetry breaking; there is an approximate 
SU(2)l x SU(2)r symmetry which is spontaneously broken to SU(2)y (isospin). This 
pattern of spontaneous symmetry breaking gives rise to three Goldstone bosons which 
are identified with the pions. In the chiral limit pions interact only through derivative 
interactions. A finite pion mass reflects the nonvanishing current quark masses in QCD. 
X-PT is a useful theory because pions interact only through vertices with a countable 
number of derivatives and/or insertions of the quark mass matrix; each power of the pion 
mass and each derivative count as one power of "small momenta". Therefore S'-matrix 
elements computed in xPT can be expressed as simultaneous expansions in powers of 
momenta and pion masses over the characteristic scale of physics that is not included 
explicitly in the effective theory: 




(1) 
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taking Q to represent a "small momentum" : a momentum p -C A x or a pion mass. Note 
that the dynamical effects of all mesons other than the pion are accounted for through 
local pionic operators. In the purest form of xPT ', the coefficients of these operators are fit 
to one experiment and then used to predict other experiments. Hence xPP is a predictive 
theory. The machinery of xPT was originally developed for Goldstone boson interactions. 
However, it is straightforward to extend the machinery to include single-baryon processes, 
with the baryons treated as heavy, static objects. This subject is reviewed in great detail 
in Ref. 0. 

Nucleon Compton scattering has been studied in xPT in Ref. @, where the following 
results for the polarizabilities were obtained to order Q 3 : 



a n = a n = — ^# — = 12.2 x 10 _4 fm 3 ; 
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^"= 768^ = 1.2xlO-fm3 (3) 

Here we have used qa = 1-26 for the axial coupling of the nucleon, and f n = 93 MeV as the 
pion decay constant. Note that the polarizabilities are predictions of xPT & t this order. 
The 0(Q 3 ) xPT predictions diverge in the chiral limit because they arise from pion loop 
effects. In less precise language, the power counting of xPT implies that polarizabilities 
are dominated by the dynamics of the long-ranged pion cloud surrounding the nucleon, 
rather than by short-range dynamics. The polarizabilities should thus provide a sensitive 
test of chiral dynamics. 

At the next order in the chiral expansion, Q 4 , there are contributions to the polar- 
izabilities from undetermined parameters which must be fixed independently |3[]. These 
counterterms account for short-range contributions to the nucleon structure. This higher- 
order calculation will be discussed further below. 

Recent experimental values for the proton polarizabilities are 0] [] 



a p + /3 P = 13.23 ± 0.86igjg x l(T 4 fm 3 , 

a p - (3 P = 10.11 ± 1.74+lil x l(T 4 fm 3 , (4) 

where the first error is a combined statistical and systematic error, and the second set 
of errors comes from the theoretical model employed. We note that the value of a p + j3 p 

1 These are the result of a model-dependent fit to data from Compton scattering on the proton at 
several angles and at energies ranging from 33 to 309 MeV. 
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extracted from this multi-energy fit is consistent with the venerable Baldin sum rule 0. 
These values are also in good agreement with the chiral perturbation theory predictions. 
The neutron polarizabilities are difficult to obtain experimentally and so the corresponding 
xPT prediction is not well tested. A dispersion sum rule can relate the sum of electric and 
magnetic polarizabilities to an integral of the deuteron photo-absorption cross section; 
with model-dependent assumptions about the size of the neutron contribution to the 
deuteron photo-absorption cross section, it has been found that 

a n + (3 n = (14.40 ± 0.66) x 10~ 4 fm 3 . (5) 

Most direct information on a n has been obtained by scattering neutrons on a heavy 
nucleus and examining the cross section as a function of energy. Currently there is much 
controversy over what result this technique gives for a n . In Ref. the value [] 

a n = (12.6 ± 1.5 ±2.0) x 10~ 4 fm 3 (6) 
was obtained, which disagrees considerably with the result of the experiment of Ref. || , 

a n = (0.6 ±5.0) x 10- 4 fm 3 . (7) 

Given this discrepancy, it is important to measure the electromagnetic neutron polariz- 
abilities by other means. Compton scattering on a nuclear target is an obvious candidate. 
Quasi-free Compton scattering by the neutron in a deuteron target was reported in Ref. 
0. With the use of a theoretical model, the polarizability 

a n = (10.7± 3 3 7 ) x 10" 4 fm 3 (8) 

was extracted. One difficulty with this extraction of the polarizability by scattering on a 
nuclear target lies in the absence of a large neutron contribution: for quasi-free scattering 
the cross section in the forward direction goes as a n ~ (a n u 2 ) 2 (where we have neglected 
the smaller magnetic polarizability) and is very small. Partly for this reason, attention has 
turned to other methods of extracting the neutron polarizabilities. It should be pointed 
out though that some authors argue that quasi-free photon-neutron scattering deserves 
renewed study [[Uj . 

One alternative to this quasi-free Compton scattering is coherent photon scattering on 
the deuteron. In this case, the cross section in the forward direction naively goes as: 



da 



~ {hh - K + «nV ) • (9) 

0=0 



2 The values measured in these experiments are static polarizabilities of the neutron, rather than the 
polarizabilities (often denoted a) we have been discussing here. In order to correct for this we have 
applied the standard 5% correction to obtain a's from the experimental results 0. 
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The sum a„ + a n may then be accessible via its interference with the dominant Thomson 



term for the proton, Jth [jrTJ . This is one example of a general argument that the deuteron 
differential cross section is sensitive to isoscalar combinations of the nucleon polarizabili- 
ties. This means that with experimental knowledge of the proton polarizabilities it may 
be possible to extract those for the neutron. Coherent Compton scattering on a deuteron 
target has been measured at = 49 and 69 MeV by the Illinois group |T2|. An experi- 
ment with tagged photons in the energy range E 1 = 84.2 — 104.5 MeV is under analysis 



at Saskatoon [[Uj, while data for of about 60 MeV is being analyzed at Lund |Ll 

There is clearly a need for a consistent theoretical framework in order to interpret this 
data. The amplitude for Compton scattering on the deuteron involves mechanisms other 
than Compton scattering on the individual constituent nucleons. Hence, extraction of 
nucleon polarizabilities requires a theoretical calculation of Compton scattering on the 
deuteron that is under control in the sense that it accounts for all mechanisms to a given 
order in a systematic expansion in a small parameter. There exist a few calculations of this 
reaction in the framework of conventional potential models [ ITS] , [IB], 0] . These calculations 
yield similar results if similar input is supplied, but typically mechanisms for nucleon 
polarizabilities and two-nucleon contributions are not treated consistently. We will see 
that xPT provides an alternative framework where this drawback can be eliminated. 

In few-nucleon systems, a complication arises in xPT due to the existence of shallow 
nuclear bound states and related infrared singularities in A-nucleon reducible Feynman 
diagrams evaluated in the static approximation [llTl . Physically, the fundamental problem 
is that nuclear physics introduces a new mass scale, the nuclear binding energy, which 
is very small compared to a typical hadronic scale. One way to overcome this difficulty 
is to adopt a modified power counting scheme in which xPT is used to calculate an 
effective potential which generally consists of all ^4-nucleon irreducible graphs. The S- 
matrix, which includes all reducible graphs as well, is then obtained through iteration by 
solving a Lippmann-Schwinger equation ||17|| . (We will refer to this version of the effective 
theory as the Weinberg formulation.) To date the Weinberg formulation can be carried 



through rigorously only using finite cutoff regularization ]Tj|, [19], £1]]. This limitation 



has spawned an intense theoretical effort geared at formulating an effective field theory 
for low-lying bound states which is verifiably (analytically) consistent in the sense of 
renormalization |22[ . One result of this effort is a new power counting scheme in which all 
nonperturbative physics responsible for the presence of low-lying bound states arises from 
the iteration of a single operator in the effective theory, while all other effects, including 
all higher dimensional operators and pion exchange, are treated perturbatively [^, |24j . 
(We will refer to this version of the effective theory as the Kaplan-Savage- Wise (KSW) 
formulation.) This is relevant here because Compton scattering on the deuteron has 
been computed to next-to- leading order in the KSW formulation [p5fl . We will discuss 
this result and its relation to our calculation. A comprehensive and up-to-date review of 
nuclear applications of effective field theories can be found in Ref. [E6j. 
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In this paper we compute Compton scattering on the deuteron for incoming photon 
energies of order 100 MeV in the Weinberg formulation. We use baryon xPT to compute 
an irreducible scattering kernel to order Q 3 , which is then sewn to external deuteron 
wavefunctions as in Fig. |l|. Compton scattering on the nucleon is a fundamental ingredient 
of our calculation. This is a direct consequence of the power counting scheme specific to the 
Weinberg formulation. It is important to realize that typical nucleon momenta inside the 
deuteron are small — on the order of V MB or m n , with B the deuteron binding energy — 
and consequently, a priori we expect no convergence problems in the xPT expansion of 
any low-momentum electromagnetic or pionic probe of the deuteron. This method has 
led to fruitful computation of the pion-deuteron scattering length ]17|, |27j, neutral pion 



photoproduction on the deuteron at threshold p8|| , as well as pn radiative capture |2£ 
and the solar burning process pp —>■ de + v |[30|| . Although in principle we could use 
wavefunctions computed in xPT , we will consider wavefunctions generated using two 
modern nucleon-nucleon potentials. In previous computations using this method results 
are not very sensitive to details of the deuteron wavefunction. Any wavefunction with the 
correct binding energy gives equivalent results to within the theoretical error expected 
from neglected higher orders in the chiral expansion. Similar wavefunction independence 
is found here. 




Figure 1: The anatomy of the calculation. The irreducible 
kernel (I) is computed in baryon xPT and sewn to external 
deuteron wavefunctions to give the matrix element for Comp- 
ton scattering on the deuteron. 

In Section ^| we review those aspects of xPT that are directly relevant to our calcu- 
lation, including the power counting scheme and the relevant operators. In Section |3| we 
present the power counting for Compton scattering on the A A system. We discuss sub- 
tleties which arise because of the presence of a new mass scale in the two-nucleon system: 
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the deuteron binding energy. Then, in Section ^ we present our calculation of Compton 
scattering on the NN system to order Q 3 in baryon xPT . We discuss Compton scatter- 
ing on a single nucleon — which has been calculated elsewhere Jl], ||, |3| — and compute the 
Feynman amplitudes which contribute to photon scattering on the two-nucleon system 
at order Q 3 in the momentum expansion. In Section [| we present results for differen- 
tial cross sections for Compton scattering on unpolarized deuteron targets. We discuss 
higher-order effects and the extent to which a determination of neutron polarizabilities is 
possible. Finally, in Section |6] we summarize our work, and identify future directions for 
the theoretical study of this reaction. 

2 Baryon Chiral Perturbation Theory 

In this section, we briefly discuss the effective chiral Lagrangian underlying our calcula- 
tions and the corresponding power counting. Precise statements about QCD dynamics 
at low energies can be made only where a small dimensionless expansion parameter is 
identified. This is of course the motivation behind the ongoing intense effort to develop 
a perturbative theory of nuclear interactions [^IJ . One class of dimensionless parameters 
consists of pure numbers, such as a coupling constant associated with a renormalizable 
interaction, or the inverse of the dimensionality of a group, as in the large- N c expansion 
in QCD. A second class of dimensionless parameters consists of ratios of dimensional 
quantities which are present in a physical system when there is a clean hierarchy of scales. 
Effective field theory is the technology which develops a hierarchy of scales into a per- 
turbative expansion of physical observables. In a system with broken symmetries this 
technology is especially powerful. When a continuous symmetry is spontaneously broken 
there are always massless Goldstone modes which dominate the low-energy dynamics. 
Although the symmetry in question is no longer a symmetry of the vacuum, it remains 
a symmetry of the Lagrangian. Moreover, massless Goldstone modes only couple deriva- 
tively. Hence at energies small relative to the characteristic symmetry breaking scale, 
the interactions of the Goldstone bosons can be ordered in an effective Lagrangian which 
is constrained by chiral symmetry and in which each operator contains a nonvanishing 
number of derivatives acting on the pion fields. Observables computed from the effective 
Lagrangian are therefore power series in momenta, with the non-analyticities required by 
perturbative unitarity. 

In QCD the chiral SU(2)l x SU(2)r symmetry is spontaneously broken and the low- 
energy effective field theory is xPT. Here we are interested in processes where the typical 
momenta of all external particles is p <C A x , so we identify our expansion parameter 
as p/A x . If the chiral symmetry of QCD were exact this would be the only expansion 
parameter in xPT . However, in QCD SU(2)l x SU(2)r is softly broken by the small 
quark masses. This explicit breaking implies that the pion has a small mass in the low- 
energy theory. Since m 7T /A x is then also a small parameter, we have a dual expansion in 
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p/A x and m n /A x . We take Q to represent either a small momentum or a pion mass. A 
generic contribution to a matrix element involving the interaction of any number of pions 
and A nucleons can then be written in the form 



M = QTiQ/fi), 



(10) 



where \i is a renormalization scale, T{x) is some non-analytic function, which is assumed 
to obey Ti\) = 0(1), and v is a counting index. In general the full Ai is the sum of a 
set of Feynman diagrams in a quantum field theory. It is straightforward to arrive at a 
general formula for v by considering the momentum-space structure of generic Feynman 



First, we consider only irreducible subdiagrams. These are graphs in which the energies 
flowing through all internal lines are of 0(Q). The issue of calculating the full matrix 
element for the process of interest from this kernel will be taken up shortly. In irreducible 
diagrams we have the following counting rules: 

• A nucleon propagator contributes Q" 1 ] 

• A pion propagator contributes Q~ 2 \ 

• Each derivative or power of the pion mass at a vertex contributes Q; 

• Each loop integral contributes Q A ; 

• Each delta function of four-momentum conservation contributes Q~ A . 

In this way, we find that a diagram with A nucleons, L loops, C separately connected 
pieces, and Vi vertices of type i gives [[[7], 



with di the sum of the number of derivatives or powers of and fi the number of nucleon 
fields. 

Since pions couple to each other and to nucleons through either derivative interactions 
or quark masses, there is a lower bound on the index of the interaction: Aj > 0. Hence 
the leading irreducible graphs are tree graphs (L = 0) with the maximum number of 
separately connected pieces (C = max), and with vertices with A, = 0. Higher-order 
graphs (with larger powers of v) are perturbative corrections to this leading effect. In this 
way, a perturbative series is generated by increasing the number of loops L, decreasing 



rules JT7 . 



u = A-A-2C + 2L + 5^Vi A 




where the so-called index of the interaction i is 



A l = d i + fi/2 - 2 



(12) 



S 



the number of connected pieces C, and increasing the number of derivatives, pion masses 
and/or nucleon fields in interactions. 

How is this analysis altered in the presence of an electromagnetic field? Photons couple 
via the electromagnetic field strength tensor and by minimal substitution. This has the 
simple effect of replacing a derivative by a factor of the charge e. This introduces a second 
expansion, in the small electromagnetic coupling a em = e 2 /4ir. Since we will be working 
at a fixed order in the expansion in a em (leading order), it is convenient to enlarge the 
definition of di to include powers of e as well, thus continuing to classify interactions 
according to the index A, defined by (|T^). 

With the power counting scheme established, the next step is to construct the effective 
Lagrangians labelled by the interaction index Aj. The technology that goes into build- 
ing an effective Lagrangian and extracting the Feynman rules is standard by now and 
presented in great detail in !]. 

The pion triplet is contained in a matrix field 



7T 2 7T ■ T 



f-Tl +l —' (13) 

where /„- = 93 MeV is the pion decay constant. Under SU{2) L x SU(2) R , E transforms as 
E -> an d f as f L^W = U£F$; here L{R) is an element of SU(2) L (SU(2) R ), and 

U is defined implicitly. It is convenient to assign the nucleon doublet N the transformation 
property iV — > UN. 

With Q = (1 + r 3 )/2 we can write the pion covariant derivative as 



D M E = 0„£-i6A M [Q,£]. (14) 



Out of £ one can construct 



V, = -\iKd,-ieA,Q)i + i{d,-ieA,Q)e\ (15) 

A, = ^(d»-ieA»Q)Z-^-ieA»Q)?\, (16) 

which transform as -> UV^W + Ud„W and -> UA^W under SU(2) L x SU(2) R . 
is used to build covariant derivatives of the nucleon, 

D li N = (d li + V fl )N. (17) 

The electromagnetic field enters not only through minimal coupling in the covariant 
derivatives, but also through 



U» = e(£ f Q£ + £Q?)F t 
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where F pv is the usual electromagnetic field strength tensor. 

Because the nucleon mass is large, M 3> Q, it plays no dynamical role: nucleons are 
non-relativistic objects in the processes we are interested in. The field N can be treated 
as a heavy field of velocity v in which on-mass-shell propagation through exp(iMv ■ x) 
has been factored out. In the rest frame of the nucleon, the velocity vector t> M = (1,0). 
The spin operator is denoted by S M , and in the nucleon rest frame S 1 * = (l/2)(0, a). 

With these ingredients, it is straightforward to construct the effective Lagrangian, 
conveniently ordered according to the index Aj of Eq. (p~2j), 

oo 

£=Y^£ {A) . (19) 

A=0 

The leading-order Lagrangian is 

£<°> = \fttr{D^D»Y) + + S f ) 

+iN\v ■ D)N + 2g A N\A ■ S)N 
- l -C a {N^Y a Nf 

where T a is an arbitrary (non-derivative) Hermitian operator, and qa 
parameters. The sub-leading terms are contained in 

£W = ^ N ^- d2 + ( v - D ) 2 + 2 *9a{v-A,S-D} 

-l[S» t + Kv )f^ u + I(« a _ K v )trf, u })N + . . . , (21) 

where k v = k p — K n and k s = k p + K n are parameters related to the anomalous magnetic 
moments of the proton, k p = 1.79, and neutron, K n = —1.91. In sub-subleading order, 

£(2) = ~Z&T^ VPaF " uFp(T + ^P Z{Z + 2K)N ^ ^ Ax ^ N + --- ( 22 ) 
where the first term is from the anomaly, with 60123 = 1, and the second term is a 1/M 2 
correction expressed in the nucleon rest frame with electric field E = 8qA [f2|. In our 
notation, k — k p when Z = 1 and k = n n when Z = and in the expressions above, 
the ellipses signify that we have included only the interactions which are relevant to our 
calculation. Likewise, terms in £( n - 3 ) do not appear in Compton scattering to the order 
that we are working. 

Note, finally, that the A isobar can be introduced in a manner analogous to the nu- 
cleon field. The delta-nucleon mass difference Ma — M is not much larger than our Q, so 
integrating out the delta may worsen the convergence of the low-energy expansion consid- 
erably. Of course, the effects of the delta still appear in the effective theory constructed 
above, via their inclusion in counterterms, but keeping an explicit delta shifts these ef- 
fects to lower orders in the expansion. For simplicity, in this initial calculation we do not 
include the delta. 



(20) 



1.26 and C a are 
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3 Compton Counting 



3.1 The role of nuclear physics 

The power counting of chiral perturbation theory summarized above works very well for 
processes involving one or zero nucleons. A non-trivial new element enters the effective 



theory when we consider systems of more than one nucleon [17], 0] . Because nucleons are 
heavy, contributions from intermediate states that differ from the initial state only in the 
kinetic energy of nucleons are enhanced by infrared quasi-divergences. This is linked to the 
existence of small energy denominators ~ Q 2 /M, which generate contributions 0(A x /Q) 
larger than would be expected from Eq. (|TT|). However, this is not really either a surprise 
or a problem. Eq. ( |Tl~D was derived by assuming that all nucleon propagators scale as 1/Q. 
This is correct for irreducible subdiagrams, which do not contain intermediate states with 
small energy denominators. In an A-nucleon system which is not subject to any external 
probes these diagrams are A-nucleon irreducible diagrams, the sum of which we call the 
potential V. (See e.g. the left-hand graph of Fig. |2|.) To obtain the full two-nucleon 
Green's function these irreducible graphs are iterated using the free two-nucleon Green's 
function, which can, of course, have a small energy denominator which does not obey 
the xPT power-counting. This iteration gives contributions to the two-nucleon Green's 
function such as the graph shown on the right-hand side of Fig. |2|. When we consider 
external probes with momenta of order Q, the sum of irreducible diagrams forms a kernel 
K for the process of interest. The full Green's function for the reaction is then found 
by multiplying this kernel K by two-particle Green's functions in which the two-particle 
propagators with small energy denominators may appear. 



Figure 2: One- and two-pion exchange graphs. One-pion 
exchange is in the potential and (uncorrelated) two-pion ex- 
change is generated by solving the Schrodinger equation. 

In the case of Compton scattering on the two-nucleon system the formal version of 
these statements is as follows. The two-particle Green's function, G, is constructed from 
the two-nucleon-irreducible interaction V, and the free two-nucleon Green's function Go, 
via 
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G = G + G VG. (23) 

Since V is two-particle irreducible we may perform a chiral power counting on it as per 
Eq. (fLip. Such a potential has been constructed up to v = 3, in Refs. |18], |19| . That 



work employed cutoff regularization, and the Schrodinger equation was solved to obtain 
the NN wave function. Fits to deuteron properties and scattering data up to laboratory 



energies of 100 MeV were obtained for three different cutoff parameters ||20|| . However, the 
quality of the fit is not as good as that of phenomenological models. Full consistency can 
be achieved by treating both the kernel for the interaction with external probes and the 
NN potential in the same framework. However, here we will follow a hybrid approach: we 
will calculate the kernel in xPT but use several "realistic" potential models to calculate 
the two-particle Green's function. 

To calculate Compton scattering on the two-nucleon system we must also identify the 
two-nucleon irreducible kernel, K^, for the process 'jNN — > ^NN. Again, the power- 
counting of Eq. ( [TT]) applies to this object, since all nucleon denominators in it are, by 
definition, of 0(Q), as long as u ~ Q. The full Green's function for Compton scattering 
on the two-nucleon system is then 

Gyy — GK-y^G . (24) 

It follows that if i^ 77 is constructed in a gauge invariant way and is worked out to the 
same order in xPT as the potential V then this calculation will be gauge invariant up 
to that order in chiral perturbation theory. However, even if K ni and V are consistently 
constructed this calculation is not exactly gauge invariant, since a brief examination of 
the counting explained above shows that contributions to i^ 77 of the form 

K^GK^ (25) 

where K 1 is the two-nucleon-irreducible kernel for one photon interacting with the two- 
nucleon system, only appear order- by-order in the expansion in Q/A x . In fact, for gauge 
invariance to be exact the contribution (EH) must be included in full in G 77 . 



Indeed, at low photon energies the contribution (25) , known as the resonance con- 
tribution to Compton scattering, is crucial to a correct understanding of the reaction. 
In particular, without it the Thomson limit amplitude for Compton scattering on the 
deuteron will not be recovered. In this low-energy region the power-counting of Eq. ( fTT|) 
does not apply to K 77 , since graphs in K 77 with two-nucleon intermediate states (Eq. 



(]25|)) can lead to violations of the power-counting via the same infrared quasi-divergences 
discussed above. However, if the energy of the probe is larger (~ m^) it is valid to look at 
the contributions K^GK^ using the power-counting of Eq. (|ll]) . Hence at higher photon 
energies it is completely correct to treat this resonance contribution only perturbatively 
in the Q expansion. In other words, as we move into the energy region which is well 
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above the poles of the Green's function G, it becomes reasonable to use a perturbative 
expansion for G in Eq. (p5|) , rather than one which retains its full structure. 

Arguments such as these lead us to identify two regimes for Compton scattering on 
nuclei, depending on the relation of the energy of the photonic probe, uj, to the typical 
nuclear binding scale B ~ m 2 /M: 

1. uj ~ B. In this region the power-counting of Eq. (JTTJ) is not valid, since the external 
probe momentum flowing through the nucleon lines is of order Q 2 /M, rather than 
order Q. It is in this region that the Compton low-energy theorems are derived. 
Therefore our power counting will not recover those low-energy theorems. 

2. uj ~ Q 3> B. In this region the power counting of Eq. (|TTD is valid, since in all 
irreducible graphs a photon energy of order Q flows through the nucleon lines. 

These two regimes were identified by Friar and Tomusiak in a discussion of the role of 



low-energy theorems for threshold pion photoproduction on nuclei [32]. They were also 



discussed in the work of Chen et al. p5| . In Ref. Compton scattering on the deuteron 



was computed to the same order discussed here, one order beyond leading non-vanishing 
order. However, there the calculation was performed in the KS W formulation of two- 
nucleon effective field theory, rather than in the Weinberg formulation. An advantage of 
KSW power-counting is that the effective field theory moves smoothly between Q < B 
and Q > B. KSW power-counting is valid for nucleon momenta Q < A^v ~ 300 MeV. 
Thus in the KSW formulation deuteron polarizabilities and Compton scattering up to 
energies uj < A 2 NN /M ~ 90 MeV can be discussed in the same framework. Here we are 
interested mostly in the region uj ~ m n , and so we regard ourselves as being firmly in 
the second regime — although we shall see that this requires us to be cautious in our 
conclusions for Compton scattering at lower energies. In Fig. |3| we illustrate the ranges 
of validity of the KSW and Weinberg formulations of effective field theory of Compton 
scattering on the deuteron. 



3.2 Power counting for jNN —> "yNN 

Let us now consider in more detail the kernel in the regime uj ~ m^. By definition, 
it consists of those irreducible subdiagrams to which the external photons are attached. 
One cannot cut these diagrams in two pieces without cutting at least one line with energy 
0(Q). In particular, there are no irreducible graphs where the two photons are attached 
to different nucleons unless the two nucleons are connected by some interaction: energy 
and momentum conservation require the exchange of a pion or a short-range interaction 
to carry the energy 0(Q) from one nucleon to the other, and somewhere along the nucleon 
lines there is an energy flow of the order of the pion mass. 

We will concentrate here on real Compton scattering on an A = 2 system and work 
in the Coulomb gauge where the incoming and outgoing photon momenta k = (uj,k) and 
k' = (uj', k') and polarization vectors e and e' satisfy: 

k 2 = 0; k' 2 = 0; v ■ e = 0; v ■ e' = 0. (26) 
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B m|/M A 2 nn/M M A - M A * 
h 1 • •i n 1 1 1> 

(0 

i 1 

KSW 



A -less Weinberg 
i 1 

A -ful Weinberg 

Figure 3: Hierarchy of scales and effective field theories of 
Compton scattering as a function of photon energy u, as dis- 
cussed in the text. The dots represent the experimental data 
from Illinois and the cross represents the forthcoming data 
from Saskatoon. 

u = -2 

If we were interested in virtual photons, then the first contributions would start at 
v = —2. The leading effect would be from tree graphs (L = 0) with the maximum 
number of separately connected pieces (C = 2) constructed solely out of interactions with 
Aj = (X^V^Aj = 0). The only possible graphs of this type come from the photon-nucleon 
interaction generated by minimal coupling in £(°\ For real photons these graphs vanish 
in the Coulomb gauge. 

u = -1 

Likewise, some contributions at v — — 1 are still tree graphs (L = 0) with the maximum 
number of separately connected pieces (C = 2), with one photon vertex of index one and 
one photon vertex of index zero (Xi^Aj = 1). Graphs of this type do not contribute to 
real Compton scattering, if the Coulomb gauge is used. Similarly, all higher orders will 
include several diagrams that vanish for real photons in Coulomb gauge. From now on 
we do not discuss such graphs. 

Thus, at this order the only contribution comes from a tree-level diagram with two 
separately connected pieces and a vertex which is the two-photon seagull of Aj = 1. This 
vertex arises from minimal coupling in the kinetic nucleon terms in C^ 1 ' of Eq. (|2l|) . It is 
depicted in Fig. f|(a). 

v = 
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(d) (e) 



Figure 4: Characteristic one-body interactions which con- 
tribute to Compton scattering on the deuteron at order Q 2 
(a) and at order Q 3 (b-e) (in the Coulomb gauge). Crossed 
graphs are not shown and not all loop topologies are shown. 
The sliced blobs are vertex insertions from The sliced 

and diced blobs are vertex insertions from C^. 

There are four types of corrections at this order: 

1. Tree-level graphs [L = 0) with two separately connected pieces (C = 2) and two 
vertices with Aj = 1 (X^VjA, = 2). Interactions with Aj = 1 arise in C^ 1 ' from 
minimal coupling in the kinetic nucleon term and from the Pauli terms associated 
with the anomalous magnetic moments. See Fig. §(b). 

Note that if these calculations are performed in the jNN center-of-mass frame then 
the boost of the 'jN interaction from the 'yN to the 'jNN cm. frame appears 
naturally in the calculation. It occurs as a modification to these s and w-channel 
pole graphs which must be included because the total 'jN system three-momentum 
is non-zero. 

2. Tree-level graphs (L = 0) with two separately connected pieces (C = 2) and one 
interaction with Aj = 2 (Ei^Aj = 2). Here a vertex in the of Eq. ( p2|) which is 
a relativistic correction to the magnetic interaction of the nucleon with the photons 
contributes. (See Fig. f§(c).) There is also a process involving exchange of a tt° with 
the 7r°77 interaction coming from the anomaly term in C^ 2 \ (See Fig. f|(d).) 

3. One-loop graphs (L = 1) with two separately connected pieces (C = 2) and all 
interactions from the C {0) of Eq. @ (Ei^Aj = 0). (See Fig. |(e).) 
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4. Tree graphs (L = 0) with only one separately connected piece (C = 1) and all 
interactions from (See Fig. ||.) 





(c) 




(d) (e) 

Figure 5: Two-body interactions which contribute to Comp- 
ton scattering on the deuteron at order Q 3 . Permutations are 
not shown. 

Note that we have not explicitly mentioned the delta anywhere in this discussion of our 
power counting. In the energy range of interest here, the delta could give important contri- 
butions. In fact, with the delta "integrated out", as has been done here, its contributions 
first appear at v — 1, through counterterms. If the delta is included explicitly [R^, |3H| , 
then it enters at order v = via the loop graphs discussed under 3 above. 



3.3 The breakdown of the power counting 

The identification of the relevant contributions can be extended to higher orders in an 
obvious way. An example of a v = 1 (Q A ) contribution is shown in Fig. || We can use 
this graph to illustrate the transition to the very- low energy regime Q ~ rn^/M. It is 
easy to see that this graph becomes comparable to the order Q 3 graph of Fig. ||(a) when 



uM 



(27) 



Here p is a typical nucleon momentum inside the deuteron and u is the photon energy. 
Since our power counting is predicated on the assumption that all momenta are of order 
m n , we find that our power-counting formula (11) is valid for the kernel in the region 



16 




Figure 6: Two-body interaction which contributes to Comp- 
ton scattering on the deuteron at order Q 4 . The sliced blob 
represents a 1/M correction vertex from 



■jf « Q « A x . (28) 

In this first paper we will limit ourselves to a v = (Q 3 ) calculation in the deltaless 
theory. We plan to study effects of an explicit delta field in a future publication. It is 
remarkable that to this order no unknown counterterms appear. All contributions to the 
kernel are fixed in terms of known pion and nucleon parameters such as m n , gA, M, and 
/„.. Thus, to this order xPT makes predictions for Compton scattering. 

4 The xPT calculation 

In this section, we outline how the various contributions are calculated. We can sep- 
arate the contributions discussed in the previous section into two classes: one-nucleon 
and two-nucleon mechanisms. One-nucleon diagrams are those which relate directly to 
contributions to Compton scattering on a single nucleon, while two-nucleon mechanisms 
contribute only in A > 2 nuclei. One-nucleon scattering diagrams are traditionally la- 
belled as Q n with n the index in the one-nucleon process. We follow this usage here and 
above to make closer contact with the one-nucleon calculations. Of course only relative 
orders are physically significant. 

4.1 One-nucleon contribution up to order Q d {y < 0) 

Let us consider first the one-nucleon terms; a sample of the relevant diagrams is shown 
in Fig. [|. In the jN cm. frame the nucleon amplitude can be written as 

T lN = e 2 {Aie" • e + A 2 e' ■ ke- k' + iA 3 a ■ (e" x e) + iA 4 a ■ (k' x k) e ' • e 
+iA b a ■ [(?' x k) e ■ k' -{ex k') e" ■ k] + iA & a ■ [(?' x k') e ■ k' - (ex k) e' ■ k}} .(29) 
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The contribution at v — — 1 (0(Q 2 )) yields the low-energy theorem for Compton scatter- 
ing from a single nucleon — the Thomson limit: 



2„2 



T- 



-e ■ e 



M 



(30) 



where Z is the nucleon charge. The tree graphs discussed under 1 and 2 for v = in 
Section [3.2| have been known for a long time. The amplitude including these two sorts 
of contributions corresponds to a non-relativistic expansion of nucleon Born terms in a 
relativistic model with derivative coupling of the photon and Pauli magnetic moment 
interactions, plus a pion pole diagram. The differential cross section for scattering on the 
nucleon that results from these diagrams in the qa — limit is known as the Powell cross 
section. Adding diagrams of type 3 completes the 0(Q 3 ) amplitude in xPT without an 
explicit delta field 0, §, 0. 

Defining T = u/m^ and t = — 2T 2 (1 — cos 8), where 9 is the center-of-mass angle 
between the incoming and outgoing photon momenta, one finds |], [| j3|: 



A, 



A, 
A 5 
A 6 
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h{r,t) 



i . r 1 . z(l-z) . Tz 



o Jo 



S R 



with 



I^iT^t) = dx dz arcsin — -, (32) 

Jo Jo S R 

S = yjl - T 2 z 2 - t{l - z) 2 x{l - xj, R = yjl-t{l- zfx{l-x). (33) 



Because there is no counterterm consistent with chiral symmetry and gauge invariance 
at 0(Q 3 ), the sum of all loop graphs at this order is finite. Note that the loops con- 
tribute the same for protons and neutrons. Note also that, partly because of these loops, 
the Compton scattering amplitude is non-analytic in u and in cos 6. For uo sufficiently 
small a power- series expansion in uj can be made. Up to an overall factor of a em the 
coefficient of the term proportional to u 2 in the invariant amplitude A\ is then a^, and 
that proportional to u; 2 cos# is (3n- In this way we can extract the xPT predictions for 
the polarizabilities which were displayed in Eq. fl3|) 0. These effects of the pion cloud 
dominate over short-range effects represented by higher-order contact terms. The lack 
of isospin dependence in loops then implies that in general the isoscalar polarizabilities 
are larger than the isovector ones. Furthermore, a numerical factor of 10 decreases the 
magnetic polarizabilities. Thus, in general we have, a p ~ a n ^> f3 p ~ (3 n . The agreement 
with the experimental values (f|) and (0) is very good. 

The difference between the full 0(Q 3 ) amplitude and the amplitude truncated at the 
level of the polarizabilities increases with energy. At E 1 ~ 100 MeV the typical error 
is 10%. However, the polarizability approximation favors forward directions, since at 
forward (backward) angles it overestimates (underestimates) the full amplitude. Two 
extreme cases are shown in Fig. |f[ Both the full amplitude and the polarizability ap- 
proximation are in relatively good agreement with data on the proton up to ~ 120 MeV 

|35| , although the agreement deteriorates with increasing energy and angle ||3^ . 



The only modification which must be made to the photon-nucleon amplitude (|29[ ) for 
use in Compton scattering on the deuteron is that it must be boosted from the jN cm. 
frame to the ^yNN cm. frame. To the order we are working this is easily accomplished, in 
one of two equivalent ways. Since only the s and w-channel nucleon-pole graphs depend 
on the nucleon momenta, these graphs can be computed in a frame in which the total 7iV 
three-momentum is non-zero (see Fig. ||), or the 0{Q 2 ) amplitude (|30|) can be boosted 
by making the appropriate substitutions to ensure that the polarization vectors e and e' 
stay orthogonal to k and k' in the new frame p6 . In both cases an additional term: 



T boost = [e • k'e'-k + 2{e- pe' -k + e-k'e'-p)] (34) 

must be added to the one-body amplitude ( [29]) when we are computing Compton scat- 
tering from the two- nucleon system at 0(Q 3 ). 
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Figure 7: Unpolarized cross section for Compton scattering 
on the proton in deltaless X-PT to 0(Q 3 ) as a function of 
the photon energy at two laboratory angles, 6l = 1°, 180°: 
polarizability approximation (dashed line) compared to full 
amplitude (solid line). 



k' p'-k'/2 -p' -k72 k' p + k/2-k' -p-k/2 




k p-k/2 -p-k/2 k p-k/2 -p-k/2 

Figure 8: Labelling of 3-momenta for 2- and 1-nucleon mech- 
anisms in the ^NN cm. frame. 
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As for higher-order mechanisms, effects of the 0(Q A ) terms in xPT on the nucleon 
polarizabilities have been studied. In Ref. || they were estimated to be considerable, 
especially in the case of the neutron magnetic polarizability: 



a p = 10.5 ±2.0 x 10 _4 fm 3 ; a n = 13.4 ± 1.5 x 10~ 4 fm 3 ; (35) 

/3 P = 3.5 ±3.6 x 10" 4 fm 3 ; (3 n = 7.8 ± 3.6 x 10~ 4 fm 3 . (36) 

Note that the uncertainties here are theoretical and arise from the use of resonance- 
saturation to estimate some of the 0{Q i ) xPT counterterms. The largest uncertainty 
concerns the effects of the delta isobar, which has a large electromagnetic transition to 
the nucleon via an Ml operator, and therefore affects (3n significantly. Consequently, the 
resonance-saturation estimate || of its tree-level contribution to a r y r )NN counterterm is 
large, but in the case of (3 P it is largely canceled by the 0{Q 4 ) tcN loop contribution. 

Attempts have been made to better determine the delta effects by inclusion of an 
explicit A field |33|, EM. This amounts to treating Ma - Mas another quantity which is 



small relative to A x . The xPT expansion is then an expansion in the ratio of p, and 
Ma — M to A x . These three small ratios are generically labelled e. This e expansion is 
certainly appropriate for the full amplitude at p ~ m n , but possibly not very efficient for 
low-energy quantities such as the polarizabilities. The result of the explicit inclusion of 
the delta field is a shift of delta effects to lower orders. To 0(e 3 ) in this modified expansion 
there are large tree and loop contributions from delta graphs to /?jv, but now there are no 
cancelations against tcN graphs. There are also considerable loop contributions to a at. 
The 0(e 3 ) result is @: 



a p = a n = 16.4 x 10" 4 fm 3 ; (37) 
P p = p n = 9.1 x 10" 4 fm 3 . (38) 

One can only hope that an 0(e 4 ) calculation will bring these values closer to the exper- 
imental ones of (f|) and (|5[). Thus far, calculations to 0(Q 4 ) in the deltaless theory and 
to 0(e 3 ) in the deltaful theory have been limited to polarizabilities and to forward am- 
plitudes, and so, as yet, in neither of these theories is there a calculation of the invariant 
functions A%—Aq for arbitrary photon scattering angle. 



4.2 Two-nucleon contributions at order Q 3 {y = 0) 

The power counting discussed in Section |3] shows that when nuclear targets are involved 
a new contribution is present at the same order [v = 0) that the polarizabilities first 
appear: the two-nucleon contributions depicted in Fig. ^ In principle, then, we cannot 
hope to extract neutron polarizabilities without a theory for the two-nucleon contribu- 
tions. Moreover, comparison of diagrams in Figs. |](e) and |5| show that the same physics 
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that produces dominant contributions to the polarizabilities generates the dominant two- 
nucleon contributions. The only difference lies in the destination of the pion emitted by 
one nucleon: if it is absorbed on the same nucleon it contributes to the polarizability; if 
it lands on another nucleon, it generates a two-nucleon contribution. In this section we 
compute the dominant two-nucleon contributions in xPT — they involve no new parame- 
ters. It is also important to notice that these 0(Q 3 ) two-nucleon Feynman diagrams are 
the same as those at 0(e 3 ) in a theory in which the A is included explicitly. Only at next 
order would we find differences. 

The two-nucleon amplitude from the diagrams in Fig. |] is 

T™ N = (t'-t 2 - T y z ) (t« + f« + &) + t<* + t (e) ) (39) 



t (a) = + (1 <-> 2) (40) 

2[cu 2 -m 2 - (p-p' + l(k + k')) 2 } 

t ® = L g • (P - P' zMz jfog • (P ~ P' + Uk-k')) +(1 „ 2) m 
2[{p-p'-l(k-k')) 2 + ml][(p-p'+^(k-k')) 2 + ml] 

£(c) = g ■ (P ~ P' ± lg g • lg • (P ~ g iM^jOj +(1^2) (42) 

[cu 2 - m 2 - (p - p> + \(k + k')) 2 ] [(p - p> + \{k - £')) 2 + <] 

tw = e • (g- g + jfeQ g • (p - r ' - f(j? - fcQK • g ^ 

[c^ 2 - m 2 - (p - p' + + £')) 2 ] [(p - V' ~ \{k - k')) 2 + ml] 

2e ■ (g - g + e' • (g - g + |g a 1 ■ (g- p' - \{k - k 1 )) a 2 ■ (g - g + |(fc - fcQ) 
[cu 2 - M 2 - (p - p' + i(£ + £')) 2 ][(p - P' ~ \$ - k')) 2 + m 2 ][(p - p' + §(£ - £')) 2 + 
+ (1 «-> 2) (44) 



(e) 



We calculate the cross section for Compton scattering on the deuteron including the 
single-scattering and two-nucleon mechanisms described above. Our calculation represents 
therefore the full v = 0, or order Q 3 , xPT predictions for Compton scattering on the 
deuteron. 



5 Results and discussion 

To calculate the amplitude for Compton scattering on the deuteron we must take the 
residue of Eq. (^) at the initial and final bound-state pole. This leads to the amplitude 

T^ d = (^ d \K^ d ), (45) 

where |^) is the deuteron wave function and we may eliminate the isospin factor by 
making use of the relation (0d|(r 1 • f 2 — T^r 2 )\ipd) = — 2. 
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In Section f|we wrote down both the one and two-body contributions to K~ n at 0(Q 3 ) 
in the chiral expansion. If we make the same separation here then we can write down the 
following integrals for T 7 ^ in the 'jd center-of-mass frame: 



+ J^^^M<p')T^ Nx Jk',k)^ M (p) (46) 

where M (M') is the initial (final) deuteron spin state, and A (A') is the initial (final) 
photon polarization state, and k (k ') the initial (final) photon three-momentum, which are 
constrained to \k\ = \k'\ — u. Note that the single-scattering 7iV T-matrix T^ c ' m ' must 
be evaluated in the 7c? center-of-mass frame. In xPT this is straightforward, as discussed 
in the previous section. Note also that all sums over the spins of the two nucleons have 
been suppressed here. 

In general, for the wave function ip we use the energy-independent Bonn OBEPQ wave 
function parameterization which is found in Ref. | 37fl . The photon-deuteron T-matrix (|46|) 



is then calculated and the laboratory differential cross section evaluated directly from it: 



da 1 (E'\ 2 \ 




dVL L 16tH \E^J 6 M'X'MX 
where E 1 is the initial photon energy in the laboratory frame, and is related to to, the 
photon energy in the 7c? center-of-mass frame, via: 

u = 5 ; (48) 

yjl + 2E 1 /M d 

and Ely is the final photon energy in the laboratory frame: 

E> = (49) 

7 M d + E y (l- 00s 6 L )' K ! 

Convergence tests indicate that with the numbers of quadratures chosen the cross section 
evaluated in this fashion is numerically accurate at about the 1% level. Of course, this 
error does not include the theoretical error from uncertainties due to different deuteron 
wave functions, and the effect of higher-order terms in the xPT calculation of the kernel. 
These errors will be discussed further below. 

The rest of this section is divided into four subsections. First we display our results 
for this straightforward 0(Q 3 ) calculation. Results are presented at laboratory photon 
energies of 49, 69, and 95 MeV. At the lower two energies we compare with the exper- 



imental data of Lucas [O]. Data at E~ = 95 MeV are expected soon from a recent 
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Saskatoon tagged photon Compton scattering experiment fl3|| . Second, we observe that 
at lower energies contributions due to NN states with relative momenta of order yMu 
lead to certain diagrams being enhanced over the value expected from the power counting 
of Section [2|. Indeed, as discussed in Section [3], we expect this power counting to break 
down for to ~ m 2 /M. We estimate the effect of these contributions by calculating the 
impact some of them have on the differential cross section for photon-deuteron scattering. 
Third, we display the sensitivity of our results to the choice of deuteron wave function. 
Re-calculating our amplitudes with ip d defined by solving the Schrodinger equation with 
the Nijm93 potential allows us to estimate the error associated with uncertainties in 
the NN interaction. Finally, in order to see the size of the effects we can expect in an 
0(Q 4 ) calculation of this process, we arbitrarily include two terms which contribute to 
the 0(Q A ) calculation: counterterms which shift the values of the nucleon polarizabilities 
from the 0(Q 3 ) result. 

5.1 Results at 0(Q 3 ) 

In figures [| ID and |TT| below we display our results at 49, 69, and 95 MeV. For comparison 
we have included the calculation at 0(Q 2 ), where the second contribution in Eq. (|46| ) is 
zero, and the 'yN T-matrix in the single-scattering contribution is given by the Thomson 
term on a single nucleon, Eq. (PU|). 

The curves show that the correction from the 0(Q 3 ) terms gets larger as lu is increased, 
as was to be expected. Indeed, while at lower energies corrections are relatively small, in 
the 95 MeV results the correction to the differential cross section from the 0(Q 3 ) terms is 
of order 50%, although the contribution of these terms to the amplitude is of roughly the 
size one would expect from the power counting: about 25%. Nevertheless, it is clear, even 
from these results, that this calculation must be performed to 0(Q A ) before conclusions 
can be drawn about polarizabilities from data at photon energies of order m w . This is 
in accord with similar convergence properties for the analogous calculation for threshold 



pion photoproduction on the deuteron [28 



We have also shown the six Illinois data points at 49 and 69 MeV ||12|| . Statistical 
and systematic errors have been added in quadrature. It is quite remarkable how well 
the 0{Q 2 ) calculation reproduces the 49 MeV data. However, it is clear that the agree- 
ment at forward angles is somewhat fortuitous, as there are significant 0(Q 3 ) corrections. 
Meanwhile, the agreement of the 0(Q 3 ) calculation with the 69 MeV data is very good, 
although only limited conclusions can be drawn, given that there are only two data points, 
each with sizeable error bars. 

The energy-dependence of our results at a given angle is due to a combination of two 
mechanisms. First, at backward angles the 0(Q 2 ) deuteron Compton scattering amplitude 
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Figure 9: Results of the 0{Q 2 ) (dotted line) and 0{Q 3 ) (solid 
line) calculations at a photon laboratory energy of 49 MeV. 
The data points of Ref. ]12] are also shown. 
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Figure 10: Results of the 0{Q 2 ) (dotted line) and 0(Q 3 ) 
(solid line) calculations at a photon laboratory energy of 69 
MeV. The data points of Ref. [O are also shown. 
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Figure 11: Results of the 0(Q 2 ) (dotted line) and 0(Q 3 ) 
(solid line) calculations at a photon laboratory energy of 95 
MeV. 



has energy- dependence, due to it is proportionality to the non-relativistic deuteron charge 
form factor at a momentum transfer 



q = k — k' . 



(50) 



At forward angles this energy-dependence is minimal, since \q\ 2 is proportional to uj 2 (1 — 
cos 8). Second, at both forward and backward angles the amplitude also has significant 
energy-dependence due to the structure of the single-nucleon amplitude fl2"9"|) as a function 
of uj. This is particularly apparent at 9 cm = 0, where the 0(Q 2 ) amplitude is, in fact, 
energy-independent, and we can see that the polarizabilities play a significant role in 
suppressing the forward differential cross section as the photon energy increases. 

Our results are qualitatively not very different from other existing calculations. At 49 
and 69 MeV our 0(Q 3 ) results are very close to those in Ref. [15| and a few nb/sr higher, 
especially at back angles, than those of Refs. |16|, [7| (which are similar at these energies). 
At 95 MeV our 0(Q 3 ) result is close to that of Ref. [|l6j, higher by several nb/sr at back 
angles than Ref. |7|, and several nb/sr lower than the calculation with no polarizabilities 
of Ref. 0] f\ Comparing to the calculations of deuteron Compton scattering in the 
KSW formulation of effective field theory J 



we see that the result of Ref. 25 



is 



significantly lower than those presented here at both 49 and 69 MeV. At 49 MeV the 



3 At this energy Ref. Jig] only presents results with a p + a n = [3 p + [3 n = 0, which in turn arc 
considerably less forward peaked than the corresponding calculation of Ref. m& . 
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agreement of Ref . pSf 's calculation with the data is better than ours. We shall show in 



the next section that this is partly because 49 MeV is at the lower end of the domain 
of applicability of the Weinberg formulation. At 69 MeV our calculation does a slightly 
better job of reproducing the (two) data points available. The qualitative agreement 
among these calculations is a reflection of the similarities of mechanisms involved. Ours 
is however the only calculation to incorporate the full single-nucleon amplitude instead 



of its polarizability approximation. As shown in Figure [12| our tendency to higher cross 
sections in the backward directions is at least in part due to this feature, as expected from 
the discussion in Sect. 13. 
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Figure 12: Unpolarized cross section for Compton scattering 
on the deuteron in deltaless xPT to 0(Q 3 ) at 95 MeV. The 
polarizability approximation to the single-nucleon amplitude 
(dashed line) is to be compared with the full 0(Q 3 ) calculation 
(solid line). 



5.2 The breakdown of the power counting at low energies 

In this section we examine the breakdown of the power counting defined by Eq. (|Tl|) at 
low energies. Recall that this power counting was only derived for photon energies of 
order m^. Here we will calculate one portion of a particular diagram, and show how 
the contribution of that diagram becomes enhanced beyond the expectations of power 
counting if the photon energy becomes too low. 

If the photon energy uj is of order m n 



Consider the diagram shown in Figure 13 



and therefore much larger than the nuclear binding and the nucleon kinetic energies then 
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the denominator of the intermediate-state NN propagator can be approximated by l/u. 
Effects due to the motion of bound nucleons then appear as higher-order vertices in the 
xPT Lagrangian, while effects due to binding also occur at higher order in the calculation. 
This is the approach adopted in this work. However, as u is decreased, replacing the 
denominator in the two-nucleon propagator by u is no longer a valid approximation. It 
becomes necessary to retain the full energy-dependence for dynamical nucleons, rather 
than use the heavy-baryon propagator which treats the nucleon as a static source. Thus, 
in this section we evaluate this graph using a propagator which retains the full momentum 
dependence for the two-nucleon intermediate state. In heavy-baryon chiral perturbation 
theory this means replacing the standard heavy baryon xPT propagator for two nucleons 
with energy uj by one in which a set of nominally higher-order corrections representing 
the motion of the nucleons has been resummed fl7| : 



u) Muj — p 2 ' 

with p the relative three-momentum of the two-nucleon state. 

To simplify matters we consider only the vertices from the xPT Lagrangian which 
arise from electric coupling of the photon to the nucleon. This gives: 




Figure 13: Contribution to Compton scattering on the 
deuteron from the s-channel pole in photon-nucleon scatter- 
ing. 



T, 



s— channel pole 



d 3 p ,(_ k-k'\ e'-(2p+k-k') M e ■ (2p) 

4>\P + — n — u ; ^2 h—^kp) 



,2 - :i V 2 / - ]f M(co-B)-(p+l) 2 + lV M 

(52) 
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To estimate the size of effects arising from the dynamical nature of the nucleons in the 
intermediate NN state we have rewritten the two-nucleon propagator as: 



M 



M 



( 



M{u - B) 



(p+l 



-iMirS 



+ ir] M(u - B) 



P + l 



Moo - MB 



\ 




(53) 

For the case of low-energy photon-deuteron scattering the leading contribution to the 
delta- function piece will be 5(Mu — p 2 ). Thus, we may evaluate the leading part of the 
delta- function piece of T s _ channe i po [ e as a two-dimensional integral over the solid angle fl p : 



Ts^channel pole = 2 ^J^J M J $ [? + g J ^ ^ ^ ( 54 ) 



Naively this expression contributes to T ld at an order \Juo / 'A x relative to the leading 
contribution. However, note that in fact this part of T s ^ c hannei pole is suppressed beyond 
the level indicated by this straightforward counting, since factors of the deuteron wave 
function at momenta of order y/Mu enter. Ultimately if u ~ these factors drastically 
reduce this piece of T s _ c h a nnei pole, rendering it quite small. However, at ui = 49 MeV, 
contributions such as (|54l) can be sizeable 0. 

A contribution to Compton scattering on the deuteron similar to that of Q53D exists for 
the w-channel pole graph for photon-nucleon scattering. The contribution of that graph to 
the amplitude is, in fact, real, and the loop momentum that contributes to that process is, 
at lowest order, \p | = iy/Mco. In fact, as explained in Refs. |25| it is the contribution from 
these momenta in the w-channel pole graph that facilitates the (approximate) recovery of 
the Thomson limit for low-energy photon scattering off the deuteron. 

When these two contributions are added to the amplitude T ld calculated above we 
obtain the results shown in Figure 0. There we have shown results for 0(Q 2 ) and 0(Q 3 ), 
both with and without these mechanisms, at laboratory photon energies of 49, 69, and 
95 MeV. Clearly the contribution from these two-nucleon states with relative momenta of 
order V Muo plays a significant role if E y = 49 MeV. Since these effects are higher-order 
in our power counting this apparently implies a breakdown of our approach at u ~ 50 
MeV. However, we observe that once this iViV intermediate-state contribution is added 
the agreement with the experimental data at 49 MeV is really quite good. This suggests 
that it may be possible to devise an alternative power counting for this regime, where an 
accounting is made not only of the way the scales u>, p and m n appear in the amplitude, 
but also of the role played by the additional scale yMu. 

4 The role of such contributions from on-shell intermediate-state nucleons in neutral pion photopro- 
duction on the deuteron was emphasized by Wilhelm |3^]. In contrast to the case of photoproduction, 
in Compton scattering the incoming photon can be arbitrarily soft, so such contributions are unarguably 
important if the photon energy is low enough. 
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Figure 14: Results of calculations at 49, 69, and 95 MeV 
both with and without the contribution from NN interme- 
diate states whose relative momenta are of order yMuj. In 
each case the dotted line is the 0(Q 2 ) calculation, and the dot- 
dashed line is the 0(Q 3 ) calculation. Meanwhile the dashed 
line is an 0(Q 2 ) calculation with the contributions from NN 
intermediate-states with momenta of 0(V ' Muo) included, and 
the solid line is an 0(Q 3 ) calculation with the same contri- 



butions included. Experimental data from Ref. jL2) are also 
shown. 
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However, this low-energy regime is not a primary concern in our work here. We have 
instead focused on u ~ m n , since it is there that the best opportunity to extract the 
neutron polarizability seems to lie. Our results show that as the photon energy is increased 
the contribution from these on-shell NN intermediate states becomes significantly less, 
because of the suppression from the deuteron wave function discussed above. By the time 
u = 95 MeV is reached this higher-order effect is small: in fact, as we shall see below, it 
is significantly smaller than other higher-order effects which could be included. 



5.3 Sensitivity to choice of deuteron wave function 

We have also performed the calculation of the photon-deuteron differential cross section 
with a wave function found by solving the Schrodinger equation with the Nijm93 poten- 
tial [ 55[ . A comparison of 0(Q 3 ) calculations with the two different potentials at 49, 69, 
and 95 MeV is shown in Figure |T3[ The results with the Nijm93 wave function are in 
general a little higher than those found with OBEPQ. However, the cross sections only 
differ at the 10% level. 

Further investigation shows that much of this discrepancy is due to inconsistency be- 
tween the physics of pion range in the NN potentials employed and the two-nucleon 
Compton scattering mechanisms involving pions which we have calculated. In computing 
these graphs we employed the ttNN vertex from the lowest-order chiral Lagrangian and 
the axial coupling constant qa = 1.26. Therefore in essence we have used a 7rNN coupling 
in these graphs given by the Goldberger-Treiman relation: 

In contrast, the values of f 2 employed in the two NN potentials from which we obtain 
our deuteron wave functions are: 

■f 2 -««L m = °- 079 ; £™L-» = °- 075 - < 56 > 

Effects beyond the leading-order result for f n NN, Eq. (|55|), will enter at orders beyond 
those we have considered here and presumably move our f n NN towards the values (|56|). 
Although such changes in f n NN would be higher order in the chiral expansion, we note that 
in certain circumstances they may enter raised to the fourth power in the photon-deuteron 
cross section, and so might play a non-negligible role. 

A fully consistent resolution of this issue must await the use of an accurate %PT 
deuteron wave function in the calculation of deuteron Compton scattering. To address the 
issue in our current calculation we multiplied the strength of our two-nucleon mechanisms 
by two different factors in the calculations employing the two different deuteron wave 
functions. In this way the two-nucleon current is made consistent with the value of the 
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ttNN coupling employed in the Nijm93 and OBEPQ potentials. We stress that the two 
potentials still do not have the same strength at pion range, but at least once such a scaling 
is employed the NN mechanisms for Compton scattering are always consistent with the 
one-pion exchange used to obtain the deuteron wave function. This partial remediation 
of the inconsistency reduced the wave function sensitivity to about 5% everywhere. This 
suggests to us that the difference between the results obtained with the two different wave 
functions does not represent a true sensitivity to short- distance dynamics, but rather a 
need to treat the dynamics of pion range in the same way in all pieces of the calculation. 

As we shall see, this wave-function dependence is smaller than other effects which are 
higher-order in the chiral expansion of the kernel. Therefore at this stage we are not 
concerned about it. Ultimately we would like to perform this calculation of Compton 
scattering on the deuteron with a deuteron wave function calculated consistently from an 
iViV potential derived within the xPT framework [2(| . 
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Figure 15: Results of calculations at 49, 69, and 95 MeV us- 
ing two different NN wave functions in the 0(Q 3 ) calculation. 
The solid line is the result using the OBEPQ wave function, 
and the dot-dashed line using the Nijm93 wave function. Ex- 
perimental data from Ref. [O are also shown. 



5.4 Effects of higher-order terms 

Finally, in order to test the sensitivity of our calculation to higher-order effects we added 
a small piece of the 0{Q A ) amplitude for Compton scattering off a single nucleon. Specif- 
ically, we modified the invariant functions A\ and A 2 of Eq. ([H]) to include additional 
terms. The additional terms are 
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AA 1 = (Aa + A(3cos9)u; 2 , (57) 
AA 2 = -A/5. (58) 

We emphasize that these are only two of a number of new terms which will appear in the 
single-nucleon scattering amplitude T 7 tv at 0(Q A ). Furthermore, a number of additional 
two-body mechanisms must be included in T™ N in any 0(Q 4 ) calculation of Compton 
scattering on the deuteron. Nevertheless, here we calculate the differential cross section 
with the terms ( 57 ) and ( pB[ ) in order to get a feel for the sensitivity of our result to the 



presence of such higher-order terms. 

The addition of (|57|) and (^) change the polarizabilities in the calculation to 

a = a m + Aa, (3 = f3 {Q3) + A/3, (59) 

where a^ Q ^ and f3^ Qi ^ are the 0(Q 3 ) values of Eq. (|3|). 

Two calculations were performed. In the first, Aa and Af3 were chosen so that the 
total polarizabilities (|59|) were equal to the "experimental" values (f|), @, and (H). The 
second calculation involved a more dramatic change in the polarizabilities: Aa and Af3 
were chosen so that a and (3 were equal to the 0{Q A ) values of Eqs. ( p5|) and (pop. In 
either case Aa p + Aa n is relatively small, while A/3 P + Af3 n is not large for "experimental" 
values, but is significant for the 0(Q 4 ) values. The results of these two calculations for 
the two photon energies E 1 = 49 MeV and E 1 = 95 MeV are shown in Fig. ^ and Fig. [17. 



In both cases we see that, just as one would expect, the cross section at 95 MeV is much 
more sensitive to these 0{Q i ) terms than the cross section at 49 MeV. It is not surprising 
that the calculation with 0(Q 4 ) polarizabilities exhibits a larger change than that with 
"experimental" values. Continually increasing (3 P + /3 n at approximately constant a p + a n 
decreases the cross section at forward angles and increases it at back angles. In fact, it 
seems that if (3 P + (3 n is sufficiently large then the character of the cross section at 95 MeV 
can change completely from forward peaked to backward peaked. 



As a final exercise, Figs. and [19] exhibit the differential cross sections found at 49 and 
95 MeV if a p and (3 P are held fixed at their 0(Q 3 ) values, while (3 n and a n are varied with 
a n + (3 n constrained to the central value of Eq. (§). These plots suggest that at backward 
angles and higher energies the differential cross section is quite sensitive to a n — f3 n . 

However, significant change in the 95 MeV results presented here from those obtained 
at 0(Q 3 ) mandates a cautious interpretation. At the same time that the cross section 
at 95 MeV is more sensitive to polarizabilities than at lower energies, it is also more 
sensitive to 0(Q A ) corrections. In our view, a full 0(Q A ) calculation in xPT is necessary 
if any attempt is to be made to extract the neutron polarizability from the Saskatoon 
data within this framework. 
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Figure 16: Results of calculations at 49 MeV using different 
values for the nucleon electromagnetic polarizabilities. The 
solid line is the result using the 0(Q 3 ) xPT value, the long- 
dashed line is the result using "experimental" polarizabili- 
ties, and the dot-dashed line represents a calculation with the 
0(Q 4 ) polarizabilities. Experimental data from Ref. 
shown. 
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Figure 17: Results of calculations at 95 MeV using different 
values for the nucleon electromagnetic polarizabilities. Leg- 
end as in Figure [IB. 
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Figure 18: Results of calculations at 49 MeV which vary the 
neutron electric and magnetic polarizabilities while holding 
the sum fixed. 
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Figure 19: Results of calculations at 95 MeV which vary the 
neutron electric and magnetic polarizabilities while holding 
the sum fixed. 
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6 Conclusion 



We have calculated the differential cross section for Compton scattering on the deuteron 
in xPT up to 0(Q 3 ). We have found: 

• Reasonable agreement with the data at 49 MeV. At this energy 0(Q 3 ) corrections 
are not large compared to the leading 0(Q 2 ) result, and 0(Q 4 ) terms seem to be even 
smaller. There is little dependence on the wavefunction used. However, as anticipated, 
certain terms necessary to restore the Thomson, zero-energy, limit for photon-deuteron 
scattering are important at this lower energy. Once they are added to our calculation 
agreement with data is good. 

• Good agreement with the data at 69 MeV. At this energy the convergence appears 
to be good. This suggests that %PT at 0(Q 3 ) is providing reasonable neutron and two- 
nucleon contributions. 

• That the polarizability approximation should not be used in the calculation of the 
differential cross section at 95 MeV, since truncating the photon-nucleon amplitude at 
order u 2 results in a significant change in the photon-deuteron differential cross section 
for forward angles. 

• Wave function dependence on the order of 10% in the differential cross section. Our 
calculations suggest that much of this dependence would be removed if the physics of pion 
range were the same in both wave functions and in the computation of the kernel. 

• A prediction at 95 MeV which is, however, plagued by considerable uncertainties. 
Convergence is slow at this energy, as indicated by the relative size of both the full set of 
0(Q 3 ) corrections and a partial set of 0(Q 4 ) corrections. The cross section tends to come 
out somewhat smaller than at lower energies, in particular in the backward directions, 
although the full 0(Q 4 ) amplitude is likely to be somewhat bigger at back angles. It 
seems that a more stringent test of xPT at these energies — including aspects of neutron 
structure beyond the 0(Q 3 ) "pion cloud" picture — will have to wait for a next-order 
calculation. 
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